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Kinks and small-scale structure on cosmic strings 
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We discuss some hitherto puzzling features of the small-scale structure of cosmic strings. We argue 
that kinks play a key role, and that an important quantity to study is their sharpness distribution. 
In particular we suggest that for very small scales the two-point correlation function of the string 
tangent vector varies linearly with the separation and not as a fractional power, as proposed by 
Polchinski and Rocha [Phys. Rev. D 74, 083504 (2006))]. However, our results are consistent with 
theirs, because the range of scales to which this linearity applies shrinks as evolution proceeds. 
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I. INTRODUCTION 
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From extensive studies over recent years a broad consensus has emerged about the evolution of a network of cosmic 
strings. The various simulations all yield basically similar results Jl, 2, 3, 4, 5, 6, 7]. The network will evolve towards 
I a scaling regime in which the characteristic length scale ^ grows in proportion to the age of the universe t. The 
energy-loss mechanism that makes this possible is the formation and slow decay of loops. However, many details, 
Q , concerning small-scale structure, the typical sizes of the loops and the mechanism by which they lose energy and 

• eventually disappear are still uncertain. We do not know for sure whether the small-scale structure will also scale. 
c/3 ' The primary length scale ^ of the network, the typical inter-string separation, may be defined by requiring that a 
randomly chosen large volume V will on average contain a length of long string (i.e. strings of length greater than 
the horizon size, excluding smaller loops) 

(N 
> 

o: c 

Here, and subsequently, 'length' of string means the invariant length — for a straight piece of string, the length in its 
rest frame. Moreover, we use physical, not comoving, length units. The assumption of scaling is that, after an initial 
transient period, ^ will increase in proportion to the time. Once scaling is achieved, the simulations tell us that 

o ■ 

■ f « 7i, with 7r « 0.3, 7,„ « 0.55, (2) 

o: 

• • ■ where the subscripts denote the radiation-dominated and matter-dominated eras. 

From these simulations, it is hard to extract reliable information about small-scale structure, because we are talking 
about scales close to their resolution limit. Moreover, the simulations can only cover a few Hubble times, while it 
is very possible that short-distance scaling will only be reached after a much longer time. For that reason, analytic 
studies are important. 

In our first attempt to construct an analytic model we introduced also a second length scale ^ defined in terms of 
the correlation function and representing the persistence length along the string, i.e., the distance beyond which there 
is no directional correlation. Later this model was extended to the 'three-scale model' of Austin et al [9'], incorporating 
a length scale ^ loosely defined as a typical separation between kinks. However, the treatment of small-scale structure 
relied on a number of untested assumptions and approximations and was not wholly satisfactory. A considerable 
advance in a similar direction has been achieved more recently in the work of Polchinski and Rocha [l^, [ll| , based 
on a study of the correlation functions on strings. They obtained very good agreement with results of simulations 
over the intermediate range of scales, but there remain puzzling features and unexplained discrepancies at very small 
scales. 
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In this note, we point out some significant issues, especially concerning the distribution of kinks on strings, and 
revisit some of the questions discussed in our earlier work [3| which we believe may have an important bearing on 
the resolution of these problems. We begin in Sec. Hi] by reviewing very briefly the dynamics of cosmic strings, and in 
Sec. lIIII the estimation of correlation functions. In Sec. II VI we discuss the calculation of the number of intercommuting 
events, and in Sec. |V] the derivation of the equation for the rate of change of the total length of long string. Then in 
Sec. IVII we study the distribution of kinks, which we argue holds the key to understanding the small-scale structure. 
The results are discussed in Sec. IVIII 



II. DYNAMICS 



For simplicity, we consider only the standard ('vanilla') cosmic strings, moving in a spatially flat Friedmann- 
Lemaitre-Robertson- Walker space-time, with metric 

ds^ = dt^ - a^{t)dx^, (3) 

where a oc , with i^r — ^ and i'^ = ^. 

The dynamics of cosmic strings on all scales much larger than the string width (and on time scales where gravita- 
tional or other radiation is insignificant) are well described by the Nambu-Goto equations. They take their simplest 
form if we use null (characteristic) world-sheet coordinates u and v, satisfying the conditions x\ = x\ — 0, where 
the subscripts denote partial derivatives. In flat space-time we can choose u,v — t± s where s measures the invariant 
length along the string (the energy divided by the string tension /x) , and even in curved space-time we have the free- 
dom to make reparametrizations (u^v) (u'(m), w'(d)), so we can always choose u and v locally in this way. When 
discussing small-scale features it is generally convenient to do so. The Nambu-Goto equations then read 

,uv — H {t ~\- t 'qX ^ (4) 

where H = a/a. The characteristics are the constant-u and constant-w curves. 

It is convenient to think of the string as a superposition of left- and right-moving components. We introduce spatial 
unit vectors representing these waves, defined by 

(XX 11 ax y . , 

P=^r^, (1 = ^—- (5) 

Then the equations of motion ([4]) imply 

Q,u = -Htu{p-pqq)- (6) 

On time scales short compared to the Hubble time, the right-hand sides are small, so these equations show that p 
varies only slowly along left-moving characteristics and q along right-moving. If we regard p as a function of u and 
the time t, and q as a function of v and i, these equations reduce to 

p = -H{q- q pp), 

q = -H{p-pqq), (7) 

where the dots denote derivatives with respect to t. In effect, the p and q vectors repel each other slightly on the unit 
sphere as they pass. 

When strings meet they normally exchange partners, or 'intercommute'. This process creates on each of the two 
strings a pair of kinks moving in opposite directions. One consequence of the stretching of strings by the universal 
expansion is a gradual diminution of the kink angle 9. At a left-moving kink the value of p changes abruptly, say 
from Pi to P2 ■ Let us define the sharpness of the kink to be 

^=5(l-Pi-P2)=sin2(0/2), (8) 

which of course ranges from to 1. Then from ([7]) we have 



ip = \H{q ■ P2 q ■ PiPi ■ P2 + Pi ■ Q - Pi ■ P2P2 ■ Q)- 



(9) 
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Since the time scale for change of 'ijj is the Hubble time, we can replace the scalar products with q by their time 
average, related to the mean square velocity by 



According to the simulations. 



Thus dSl) becomes 



whence it follows at once that 



where 



(p.g) = -a = -(1-2(^2)). (10) 



ar ~ 0.18, an, ~ 0.3. (11) 



= -2Hai^, (12) 



tP cx cx t-"^^, (13) 



C = i^a, or Cr ~ 0.09, Cm ~ 0.2. (14) 



Note that in the radiation-dominated era this blunting of kinks is an exceedingly slow process: for the sharpness to 
reach half its value at i = 1 we must wait until t — 47 (though in the matter era, t = 5.7 would suffice) . 



III. CORRELATION FUNCTIONS 

Polchinski and Rocha [lo| used these same equations to derive expressions for the p-p (or q-q) correlation function 
on the strings. In particular, they argued that the correlation function should have the form 

z{s, t) = {p{ui,t) ■ p{u2, t)) = l- A{s/tf^, (15) 

with s = wi — U2, where we have chosen the null coordinates locally to he u,v = t ± s. Here A is a constant and 

X=Y^, or Xr«0.10, Xm«0.25. (16) 

This expression provides a very good fit to the data from simulations over a considerable range of scales, but at very 
small scales, particularly in the radiation era, the correlation function seems to approach unity much faster than this 
would suggest. In fact, we shall argue that cannot be correct at very small scales. Nevertheless, as we shall see, 
it may still be true in some limiting sense. 

At the end of the friction-dominated era, the strings are expected to be smooth, with no sharp kinks. In the 
subsequent evolution, they will always remain smooth except for the kinks introduced by intercommuting events. But 
in a finite time there can only be a finite number of such events — later we shall estimate the number — so the 
strings should always be piecewise smooth. This means that on any short stretch of string of length s the probability 
of finding a kink goes to zero in proportion to s. It follows that for very small s, 1 — z{s,t) cx s; it cannot behave 
like any fractional power of s. The kink distribution is, we contend, the key to understanding the very small scale 
properties of the string network. 

As discussed in ref. [^] the equations ([7]) can also be used to relate the correlation function z to a (though some 
of the factors in that earlier discussion were incorrect). The vectors p{u,t) and q{v,t) meet when 2t = u + v — 2to, 
say. The only reason they are then correlated, albeit weakly, is of course the 'repulsion' effect. As it approaches p our 
chosen vector q is 'repelled' by other p{u') vectors, and although these small changes fluctuate in direction, so long 
as the general direction is the same, there will be a net effect. We have 

dtiP ■ q) ~ -H[{p p-p qq p) + {q q-q pp q)], (17) 

where (see Fig. [1]) p and q stand for p{u,t) and q{v,t) while 



p = p{2t — V, t), q = q{2t ~ u, t). 



(18) 
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FIG. 1: Vectors involved in the rate of change of p ■ q. 



Since the correlation between the p and q vectors is very weak, the ensemble average of the expression in the first 

3^ 



parenthesis is to a good approximation |(p • p), and similarly for the second parenthesis. Thus we find 



dt{p{u,t)-qiv,t)) ^-^H{t)z{u + v-2t,t), (19) 
which must then be integrated up to the time to from an early time ti before which the correlation was negligble: 

{p{u, t) ■ qiv, t))^~^ f dt' Hit')z{2to - 2t', t'). (20) 
■J Jti 

The largest contribution comes from the region where 2tQ — 2t' <C t', so to a first approximation wc can neglect the 
t' dependence of H and the explicit t' dependence of z, obtaining 

a = -(p-q) «^i/(t)e(t), (21) 

where ^ is the second length scale introduced in ref. defined by 

^(t) = / z{s,t)ds. (22) 

Note that the upper limit has been extended to infinity on the assumption that z vanishes for s > 2{to — ti). This 
^ represents the persistence length along the (left-moving) string — the distance it is expected to go on traveling in 
the same direction. 

This second length should also scale: 

e = 7t, (23) 

and we then have the relation 

a « y 7. (24) 

This would yield 

7r « 0.5, 7„, « 0.7. (25) 

Interestingly, this suggests that in both cases, ^ is slightly larger than ^. Note that although ^ is similar to the 
correlation length used among others by ref. |J| (there called £,) it is not identical. 



IV. INTERCOMMUTING 



Next, we review the calculation of the rate of intercommuting. We consider a large volume V and ask how many 
such events occur within a short time interval dt. The total length of string within V is L, given by ([T]). Within this 
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^dudv. 



Now let us randomly 

- dvi 



- duj and v 



thin slice, we can choose the world-sheet coordinates to he u,v = t ± s, so that dtds 
select two small sections of the string world-sheet, say with coordinates in the ranges Uj to Uj 
{j = 1,2). For the moment, we ignore the special case of direct small-loop creation, and assume that these sections 
are uncorrelated. The probability that the two will intersect is then the four-volume dil that they span divided by 
the total four-volume Vdt. Evidently, 



dn 



We can then extract factors of t^. 
([5]), obtaining 



where A is a 4 x 4 determinant. 



-9^\p.:^p^l,u^l,v^2,u^2,v\ duidvidU2dv2 



(26) 



i and rewrite this determinant in terms of the p and q vectors defined in 



A = - 

4 



dil = A.—duidvi.—du2dv2, 



1111 

Pi Ql P2 92 



(27) 



I (Pi X P2) • (91 - 92) + (Pi - P2) • (9i X Q2) I 



(28) 



The probability dpi that the first section of the world sheet will experience an intercommuting with any other 
section is obtained by integrating over all U2,V2, thus replacing ^du2dv2 by Ldt: 



A 1 

dpi = -^-duidvi, 



(29) 



in which we have used eq. ([T]) , and where A is the average value of A. Thus finally the total number of intercommmuting 
events in the volume V within the time interval dt is 



dNir 



AVdt 
2^4 • 



(30) 



The factor of 2 here is to compensate for the double counting of events. 

The average value of A was computed in ref. eq. (4.30) (there called x): assuming that the unit vectors are 
independently uniformly distributed over the unit sphere: 



A 



(Ppi (Pqi d^P2 d^Q2 ^ 
47r An Air At: 



where each d^p signifies an integral over the unit sphere. The result is 



A 



27r 
35 



0.18. 



(31) 



(32) 



We can allow for the small p-q anticorrelation in ^TU\\ . by assuming a linear probability distribution for p ■ q and 
inserting in the integrand of ([3T|) the two factors (1 — 3apj ■ qj)- This increases A marginally to 



^ 2tt ( 2a 0? 

A = — 1-1 

35 V 3 11 



which would yield 



A,. « 0.20, 



0.21. 



(33) 



(34) 



One interesting feature of A, which will be important later, should be noted. It vanishes when any pair of the four 
unit vectors Pi,P2i ^ii 92 equal, but not when they are equal and opposite. This has an important implication 
for kink formation. Because the probability of intercommuting is zero when p^ = but not when p^ = ~P2-< the 
angle distribution of newly formed kinks (computed below) is not symmetrical but is actually skewed towards sharper 
kinks. It appears that A attains its maximum value, 4/(3-^/3) — 0.77, when the four unit vectors are aligned towards 
the vertices of a regular tetrahedron, with equal angles, arccos(— -j), between each pair. 

We have so far ignored the direct formation of very small loops, which is really a separate type of process. In that 
case, there is a high degree of alignment between all the p and q vectors, so the value of A is always small, nowhere 
near its average value. Then the kinks produced are normally small-angle. 
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V. RATE OF CHANGE OF LENGTH 




Let us now consider the rate of change of the total length L of long string in the comoving region of volume V . It 
will change for two main reasons: stretching by the universal expansion, and excision of portions of string by loop 
formation. 

The rate of change due to stretching is a simple function of the r.m.s. velocity and the expansion rate p^ . namely 

(35) 

\ / stretching 

where a is given by (jlOp . 

Next we ask how many of the intercommniutings will generate loops. We expect that in the scaling regime a 
small, but fixed, fraction of the intercommutings will yield loops of substantial size. Some of those will subsequently 
reconnect to the long-string network. Many others will self-intersect and fragment, leading after a few oscillation 
periods to a final population of much smaller loops that oscillate and very gradually shrink. 

Each of the initial loops will remove a length of string that is on average some fraction of ^. Hence the total length 
removed within a time interval dt by loops that do not reconnect will be dL = —77^/^^ for some fixed rj, which must 
be determined from simulations. (This parameter 77 is related to the parameter qc introduced in eq. (4.31) of ref. [1].) 

This then gives 

? = (36) 

loop form 



Adding these two contributions the total rate of change of L will be 

(37) 

Using ll]) this can also be written in terms of f or 7 ^/t, as 



L t -ft 



2tj ^j]- {2-3iy + Ch. (38) 

This confirms that 7 will approach a constant scaling value, at which the right hand side vanishes. Using the values 
from the simulations this gives 

7]r « 0.18, r;,„ « 0.1. (39) 



VI. KINK DISTRIBUTION 



Here we consider the number of kinks within the volume V and their sharpness distribution. Each intercommuting 
creates two new left-moving and two right-moving kinks. We expect the numbers to remain equal on average, though 
excised loops may by chance contain more of one than the other. So let us consider the total number N{ijj,t)dilj of 
left-moving kinks with sharpness between ip and "0 + dip. 

Kinks are created by intercommuting, with some initial sharpness distribution g{ip). If we ignore kinks created 
when small loops are formed (for which the values of ^-nd P2 are strongly correlated), then this distribution can 
be found by averaging the function A of ([^5]) over all directions of and (This can most easily be done by 
transforming to the variables |(Pi ±^2) ^^'^ ^ili ^^2)1 ^^'^ choosing the z axis in the direction of |(Pi ~ P2); 
hiQi ^ Q2) the xz plane.) For simplicity, we ignore here the effect of the small anti-correlation of p and q, and 
assume that the q vectors are uniformly distributed on the unit sphere. The result is 

gW^—Viiil5-6iP-iP^). (40) 
This function is plotted in Fig. [21 The mean value of ip is easily found to be 

= j\gWdip^^. (41) 




FIG. 2: The initial kink sharpness distribution 



As expected from our earlier discussion, this is slightly larger than i. Of course, allowing for the kinks produced 
when small loops are formed would probably lower the average. 

Once formed, the kinks weaken, according to (|13p . Some of them are of course removed by being incorporated into 
loops. A reasonable first guess might be that the fraction of kinks removed is approximately the same as the fraction 
of total length excised. However, there may be some tendency for the regions incorporated into loops to have more 
kinks than others, so to allow for that we shall assume that in place of (|36|) we have 



loop form 

where fj is another parameter, possibly somewhat larger than rj. This should again be determined from simulations. 

Putting these three effects — kink formation by intercommuting events, blunting due to stretching, and removal 
by loop formation — together we then get 

^=^.W + f^(V^iV)-liV. (43) 

To solve the equation we also need an initial condition. In the radiation-dominated era, the proper condition is that 
at the end of friction domination there are no kinks at all. Let us suppose that that occurs at a time t^, . It is useful 
to change variable from ^ to tp* ~ {t/t*)'^'^'4'i which is the value would have had at t*. This yields the equation 

i^(V'*,0+(^-2C)Af(V*,i) = ^5| (^) (44) 

where A'^ on the left hand side now denotes the time derivative at constant Note that V'* may be larger than 
unity. In that case, we need the additional boundary condition that g{ip) — for ^ > \. Then it is straightforward to 
solve the equation and return to V' rather than Assuming that scaling has already been reached, so that 7 = 0, 
the result is 




N{^P,t) A /■* dt' 



where 




(45) 



/3 = 3 - 3z^ - - + 2C. 
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(46) 
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If we ignore the possibility that loops have more kinks per unit length than average, and set ?? = 77, this would imply 

A- « 1.1, /3,„ « 1.2, (47) 

so we may perhaps expect that the true values are a little smaller than these. Thus we obtain two separate expressions 
for N depending on whether ip is larger or smaller than 

If we again ignore the kinks generated by small-loop formation, and use (j40[) for the initial kink sharpness distri- 
bution, then it is easy to compute the distribution at a later time. It is convenient to write 5 (■)/') as 



(48) 



where k runs over the values k 



i, |, |, and the gt can be read off from We then obtain 



V 
V 



E 

k 

E 



(/3 + 2fcc)r 

if3 + 2fcC)7^ 



0+2kC 



(49) 
(50) 



For large angles, this is a scaling distribution — the right hand side of is time-independent. But the small-angle 
distribution does not scale. On the other hand, the limiting value tpc — (i*/i)^^ above which scaling holds decreases, 
albeit slowly, with time. 

These expressions are necessarily limited to the radiation-dominated era. In the matter-dominated era, we should 
use instead the sharpness distribution at the end of the radiation era as initial condition. We shall not pursue that 
calculation here, but for observational predictions it may be necessary to do so. However, it should be noted that by 
the onset of the matter era, the limiting value i/'c is already extremely small. 



VII. DISCUSSION 



The results obtained have some important implications both for observational predictions and for comparisons with 
simulations. 

At times of any interest to observational predictions, the critical value V'c is exceedingly small, so over almost the 
entire range of ip we have the solution (j49p . This is confirmation that the kink distribution over most of its range 
does reach a scaling regime. On the other hand, when it comes to comparisons with simulations, will still be quite 
large, so it is clear that in this particular respect we cannot expect the results to represent the full scaling regime. 

Note also that the number of kinks increases towards '0 = 0. This is because of the large number of kinks produced 
at much earlier times that have gradually become less sharp. The total number of kinks at any time is of course finite, 
though if (|49p were valid down to ijj — it would not be so. In fact the number of kinks per unit length, in units of t 
is 



tiVkink t f A 



/3-2C 



(51) 



which increases quite rapidly with time. If we use the values of the parameters estimated from simulations, we find 
in the radiation era 



kink 



L 



2.4 



(52) 



Note, however, that the large number of kinks mostly have very small sharpness. Taking kinks as a whole, the 
mean inter-kink distance i/iVkink does not scale. However, if we choose some fixed lower limit to the sharpness, the 
average distance between the kinks above that limit does eventually scale. The density of kinks grows with time, but 
the excess is entirely in very small-angle kinks. 

We can also now return to the question of the correlation function 1 — z{s) = 1 — (p(0) • p{s)) on the strings. On 
very small scales, the only contribution comes from kinks. A kink of sharpness ip between and s will contribute an 
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amount 2'ip, and if s is small enough that no more than one kink is likely to be found on the segment, then one simply 
has to integrate over the kink probability distribution, obtaining 



Evaluating the integral explicitly gives a very similar expression to (jSip . namely 



1 - z(s,t) 



2s ipA 



/3-4C 



(53) 



(54) 



where tp is the average sharpness of newly formed kinks, given by (|4T|) . namely |. Again using the earlier estimates 
of the parameters yields 



1 - z{s,t) « 3.3- 



By comparison the power-law expression psp |10| , gives 



0.7 



(55) 



l-zpR(s,i) «0.6(- 



(56) 



As expected, the expression (j55p is proportional to s, not to a fractional power of s. However it is not simply 
proportional to s/t, so it does not scale: the factor in brackets increases with t, and would become infinite in the limit 
t oo. This calculation makes sense only for small enough values of s/t that there is unlikely to be more than one 
kink present. As the number density of kinks, in terms of the scaling variable s/t, increases, the range of values of s/t 
decreases. To extend the calculation directly to larger values, we would need to know not just the density of kinks, 
but the extent of the correlation between the orientations of nearby kinks. 

Thus the picture we are suggesting is that at any time, the correlation function has a linear behavior for very small 
s/t, and that as a function of s/t it becomes steeper as time goes on. At the same time the range of values oi s/t over 
which this linear behavior is expected decreases. For larger values of s/t we expect that the power-law behavior of 
([56|) will apply. It is interesting to note that the point Sc at which the the two curves ((55|l and (|56p cross, and which 
may be taken as a rough estimate of where one type of behavior changes to the other, is for large times given by 



t V t 



0.9 



(57) 



In physical units, this limiting value s apparently increases with time, but only very slowly, like t'^'^ (though it should 
be noted that the error in this exponent may well be of order 0.1). 

These predictions can be compared with the simulations of 4] and, as shown in Fig. [31 they fit very well. The 
figure shows a log-log plot of 1 — z(s, t) as a function of s/t for eight different times, with values of t increasing from 
right to left. All the curves follow the Polchinski-Rocha form roughly from s/t ^ 1 down to a limiting value that 
moves down as t increases. For smaller values of s/t the curves steepen fairly sharply and at least for the later times 
fit well with the linear form (|55p . and as expected the turnover point moves leftwards with increasing time. Naturally, 
the transition between the two regimes is not sharp, but nevertheless very clear. Given that the simulations start 
from an initial state in which the strings are smooth, it is appropriate to set equal to the initial time. Then the 
latest time, represented by the leftmost (red) curve is at t/t* = 10.0. Of course the coefficient of s/t in (|55)l at the 
initial time vanishes, so there we do not expect linear behavior. For the next one or two times, the curves are still 
somewhat steeper than linear, but thereafter (|55p gives a good fit. In checking the rate at which the turnover value 
of s/t decreases, it should be noted that for the earfier times we cannot use fST)) . because the term 1 in cannot 
be neglected. At least for the later curves, the rate appears entirely consistent. 

At times of observational interest, we believe the Polchinski-Rocha formulae of ref. [l3| should hold on most relevant 
scales, though it is conceivable that the different behavior at very small scales might be relevant in some cases. It 
seems plausible that the reason for the discrepancy at small s between their results and the simulations is that while 
the simulations correctly exhibit scaling on large and indeed intermediate scales, they have not yet reached the point 
of complete scaling at very small scales. It should, however, be possible to check from simulations whether this 
prediction is correct. 

It should be noted that one thing we have not considered at all is the effect of gravitational or other radiation. This 
could change the picture, but only by making the strings smoother, not less smooth, so it could not affect our main 
conclusions. 
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FIG. 3: Values of the correlation function 1 — z{s,t) as a function of s/t. Successive curves from right to left are at times 
t/t, = 1 (dashed), 1.77 (black), 2.76 (blue), 3.96 (cyan), 5.43 (green), 7.08 (yellow), 8.94 (magenta), 10.0 (red). 
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